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21 Abstract 

22 We examine satellite-based measurements of chlorophyll solar-induced 

23 fluorescence (SIF) over the region impacted by the Russian drought and heat wave of 
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24 2010. Like the popular Normalized Difference Vegetation Index (NDVI) that has been 

25 used for decades to measure photosynthetic capacity, SIF measurements are sensitive to 

26 the fraction of absorbed photosynthetically-active radiation (fPAR). However, in 

27 addition, SIF is sensitive to the fluorescence yield that is related to the photosynthetic 

28 yield. Both SIF and NDVI from satellite data show drought-related declines early in the 

29 growing season in 2010 as compared to other years between 2007 and 2013 for areas 

30 dominated by crops and grasslands. This suggests an early manifestation of the dry 

31 conditions on fPAR. We also simulated SIF using a global land surface model driven by 

32 observation-based meteorological fields. The model provides a reasonable simulation of 

33 the drought and heat impacts on SIF in terms of the timing and spatial extents of 

34 anomalies, but there are some differences between modeled and observed SIF. The model 

35 may potentially be improved through data assimilation or parameter estimation using 

36 satellite observations of SIF (as well as NDVI). The model simulations also offer the 

37 opportunity to examine separately the different components of the SIF signal and 

38 relationships with Gross Primary Productivity (GPP). 

39 1. Introduction 

40 For over 30 years, the primary tool for monitoring vegetation globally from space 

41 has been reflectance measurements at visible and near-infrared wavelengths ( e.g ., 

42 Tucker, 1979; Myneni et al., 1997). Since 1981, there is a continuous record of the 

43 Nonnalized Difference Vegetation Index (NDVI) from the Advanced Very High 

44 Resolution Radiometer (AVHRR) series of instruments on meteorological satellites 

45 (Tucker et al., 2005). The NDVI and similar indices utilize visible and near-infrared 

46 reflectances on both sides of the so-called red-edge (their difference normalized by their 
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47 sum) and are sensitive to the amount of green biomass within a satellite pixel. These 

48 indices and related parameters have been widely used to examine spatial and inter-annual 

49 variations in vegetation and for many other applications including estimation of gross 

50 primary productivity (GPP) ( e.g ., Tucker & Sellers, 1986; Randerson et ah, 1997; 

51 Running et al., 2004). 

52 Satellite measurement of solar-induced fluorescence (SIF) from chlorophyll has 

53 emerged over the last few years as a different method to monitor vegetation globally from 

54 space {e.g., Guanter et al., 2007, 2012; Joiner et al., 2011, 2012; Frankenberg et al., 

55 2011). SIF measurements are based on the fact that a small fraction of the energy 

56 absorbed by vegetation (of the order of a percent) is emitted as fluorescence. The 

57 fluorescent emission has two peaks near 685 and 740 mn, known as the red and far-red 

58 emission features. All of the satellite measurements reported thus far have been in the far- 

59 red spectral region, where reabsorption of the fluorescence within the leaves and canopy 

60 is relatively small. 

61 Relationships between SIF, NDVI, GPP and other parameters can be understood 

62 within the context of the light-use efficiency (LUE) model (Monteith, 1972), i.e., 

63 GPP = LUE * fPAR * PAR = LUE * APAR, (1) 

64 where fPAR is the fraction of absorbed Photosynthetically- Active Radiation, and 

65 APAR=fPAR*PAR is the total amount of absorbed PAR. The amount of SIF at the top- 

66 of-canopy can be approximated in a similar fonn, i.e., 

67 SIF = 0 f * fPAR * PAR * Q c = 0 f * APAR * G c , (2) 

68 where 0f is the fluorescence yield at the membrane scale, and Q c is a radiative transfer 

69 function linking the escape of fluorescence from the top of canopy to the emission of 
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70 fluorescence at the scale of the chloroplast membranes. It is reasonable to assume that Q c 

71 remains fairly constant for repeat observations of a vegetated area made from a satellite 

72 over a limited period of time when vegetation structure is not changing. 

73 The NDVI is an indicator of potential photosynthesis or photosynthetic capacity as 

74 it is a measure of chlorophyll abundance and energy absorption that varies with abiotic 

75 conditions (Myneni et al., 1995). SIF responds linearly to changes in APAR, but this will 

76 be convolved with changes in 0f that may also be related to stress. NDVI also responds 

77 to stress by a reduction of energy absorption, and this occurs on the order of a few days 

78 (Tucker et al., 1981). 

79 If F2 C is assumed constant, and the ratio LUE to 0f also remains constant, then it can 

80 be seen from Eqs. (1) and (2) that SIF will be linearly related to GPP. Theory and 

81 measurements suggest that under strong illumination, such as natural illumination present 

82 during daytime satellite overpasses, the ratio of LUE to 0f remains relatively constant, at 

83 least for fluorescence from photosystem II (e.g., Berry et al, 2013; Porcar-Castell et al., 

84 2014). Previous studies have focused on relationships between GPP estimated from flux 

85 tower measurements and satellite -based SIF in terms of both in terms of magnitude 

86 (Guanter et al., 2014) and seasonal variations (Joiner et al., 2014). These studies have 

87 demonstrated that on a weekly to monthly time-scale, there is a high correlation between 

88 GPP and SIF. 

89 Other studies have examined relationships between remotely-sensed SIF and LUE 

90 including stress. These studies have utilized ground-based measurements {e.g., Louis et 

91 al., 2005; Meroni et al., 2008; Middleton et al., 2009, 2011; Damin et al., 2010; Daumard 

92 et al., 2010) as well as satellite-based SIF (e.g., Lee et al., 2013; Parazoo et al., 2013; 
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93 Zhang et al., 2014). The latter studies with satellite data have focused primarily on the 

94 Amazonia basin and maize and soybean croplands in the midwest US. Some of these 

95 studies show that stress, including heat and moisture stress, can manifest itself earlier or 

96 be more pronounced in SIF as compared with vegetation indices ( e.g ., Daumard et ah, 

97 2010). This can occur when there is a decrease in the 0f component of SIF rather than, or 

98 in addition to, a decrease in fPAR that would be reflected in both SIF and NDVI. 

99 In this work, we examine the relative importance of 0f and fPAR to the SIF signal 

100 in a situation of high stress: the regional drought and heat wave that occurred in western 

101 Russia due to a persistent blocking ridge over central Europe during the months June 

102 through August 2010 (e.g., Gramm, 2011). Societal impacts of this event included 

103 massive peat and forest fires, a decrease in wheat production of 20-30% relative to 2009, 

104 and an increase in death rates in nearby cities including Moscow. Because this drought 

105 and heat wave occurred over an extensive region, we can examine its effects on SIF and 

106 NDVI over areas covered with predominantly different vegetation types. This allows for 

107 an assessment of whether certain vegetation types are more or less prone to stress and 

108 damage and whether stress is observed earlier in the SIF data for different vegetation 

109 types. 

110 In addition to examining satellite data, we simulate SIF and other parameters 

111 using a global land surface model forced by observation-based meteorological fields. 

112 Within this simulation, we are able to examine the effects of the drought and heat wave 

113 on fPAR and photosynthesis. This provides further insight into the relative effects of the 

114 drought on LUE, 0f, PAR, and fPAR and demonstrates the skill of the model in 

115 predicting drought-induced anomalies. To our knowledge, this region has not yet been 
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116 examined in detail in the literature with respect to satellite-based SIF observations. 

117 

118 2. Data and methods 

119 We examine data within six regions of size 2° longitude by 1° latitude over 

120 western Russia in areas impacted by the drought and heat wave in 2010. Because the SIF 

121 signal has a lower signal to noise ratio as compared with the NDVI, we need to compute 

122 averages over spatial domains approximately this size. The individual regions were 

123 chosen because they contain various fractions of different vegetation types as shown in 

124 Figure 1. The location of each box and dominant International Geosphere Biosphere 

125 Programme (IGBP) vegetation type from the MODIS Land Cover Type Climate 

126 Modeling Grid (CMG) product for 2010 are listed in Table 1 (Friedl et al., 2010). We 

127 compute 8-day averages of various meteorological and satellite vegetation parameters 

128 throughout the growing season separately for 2010 (the drought year) and for all other 

129 years with available satellite GOME-2 SIF data (2007 to 2013 excluding 2010, hereafter 

130 referred to as the climatology). 

131 

132 2.1 GOME-2 SIF 

133 The approach to retrieve the SIF signal from space was first demonstrated by 

134 observing the filling-in of the strong oxygen A-band absorption feature (Guanter et al., 

135 2007). As this approach is difficult to implement globally, subsequent satellite retrievals 

136 utilized the filling-in of solar Fraunhofer lines surrounding the oxygen A-band (near 758 

137 and 770 nm) using high spectral resolution measurements from a Fourier transform 

138 spectrometer on the Japanese Greenhouse gases Observing SATellite (GOSAT) (Joiner et 

139 al., 2011, 2012; Frankenberg et al., 2011; Guanter et al., 2012). Later it was shown that 
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140 SIF could be retrieved at 866 nm using hyperspectral measurements from the SCanning 

141 Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on 

142 board the European Space Agency’s ENVIronmental SATellite (ENVISAT) (Joiner et 

143 al., 2012) and near 740 nm with the Global Ozone Monitoring Instrument 2 (GOME-2) 

144 on MetOp satellites (Joiner et al., 2013, 2014). While spatial and temporal variations in 

145 SIF from GOSAT and GOME-2 are comparable, GOME-2 SIF has better temporal and 

146 spatial coverage than GOSAT owing to greater sampling. We therefore use GOME-2 SIF 

147 exclusively for this study. The MetOp satellites, like ENVISAT and GOSAT, are in sun- 

148 synchronous orbits. The MetOp local overpass times are -09:30. 

149 GOME -2 is a grating spectrometer that measures backscattered sunlight in a 

150 scanning nadir-viewing geometry at wavelengths between 270 and 800 nm (Munro et al., 

151 2006). GOME-2 instruments have been launched on the European Meteorological 

152 Satellites (EUMETSAT) MetOp A and B platfonns on 19 October 2006 and 17 

153 September 2012, respectively. Here, we use data from MetOp A covering the period 

154 2007-2013. The nominal ground pixel lengths near nadir are approximately 40 km and 80 

155 km in the along- and across-track directions, respectively, with a swath of width 1920 

156 km. GOME -2 achieves global coverage in this configuration within about 1.5 days. Since 

157 15 July 2013, the GOME -2 instruments onboard MetOp A and B operate in a tandem 

158 mode. In this mode, GOME-2 onboard MetOp B makes measurements with the nominal 

159 swath width and pixel size, while GOME-2 onboard MetOp A measures in a reduced 

160 swath of 960 km and pixel size of -40 km by 40 km. 

161 GOME -2 SIF retrievals are derived for a particular viewing geometry in radiance 

162 units (mW/m 2 /nm/sr) from the filling-in of solar Fraunhofer lines in the vicinity of the 
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163 740 nm far-red chlorophyll fluorescence emission peak similar to Joiner et al. (2013, 

164 2014). The retrieval uses a principal component analysis approach with a simplified 

165 radiative transfer model to estimate atmospheric absorption, surface reflectance (varying 

166 with wavelength), and fluorescence emission. We have made several adjustments in the 

167 version 2.6 (v2.6) data set used here as compared with the approaches described Joiner et 

168 al. (2013, 2014); this reduces small biases that were present in previous versions. In v2.6 

169 we use a reduced spectral fitting window between 734 and 758 nm with a single set of 

170 principal components (PCs) derived from cloudy data over ocean, desert, and ice/snow 

171 cover to estimate the spectral structure of atmospheric water vapor absorption and 

172 instrumental artifacts. We correct for drift in the absolute instrument calibration using 

173 GOME-2 solar spectra. Finally, we apply an a posteriori correction for small biases 

174 caused presumably by straylight and dark current as discussed in Kohler et al. (2014) 

175 using data over ocean. The GOME-2 v2.6 SIF data are publicly available from 

176 http://avdc.gsfc.nasa.gov. 

177 We use v2.6 level 2 SIF retrievals in this study (pixel data as opposed to level 3 

178 gridded data sets). For the time-series analysis, we average the GOME-2 data over a 

179 particular area in 8-day bins. Uncertainties are estimated in each 8-day bin as the root 

180 sum square of the standard error of the mean. A nominal constant error of 0.15 

181 mW/rrr/nm/sr was used to account for additional errors following Joiner et al. (2014). 

182 Unlike the NDVI, SIF is sensitive to the amount of solar irradiance at the surface 

183 (equation 2). When comparing directly with NDVI, we therefore normalize SIF by cosine 

184 of solar zenith angle, a proxy for the seasonal cycle of potential surface solar irradiance, 

185 and for the Sun-Earth distance. 

186 
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187 2.2 MODISNDVI 

188 We examine three different NDVI data sets from the MODerate-resolution 

189 Imaging Spectroradiometer (MODIS) on the NASA Earth Observing System (EOS) 

190 Aqua satellite: 1) the standard MYD13Q1 vegetation indices data set (Huete et al., 2002); 

191 2) the Global Inventory Modeling and Mapping Studies GIMMS NDVI data applied to 

192 Aqua MODIS (Tucker et al., 2005); 3) MODIS NDVI computed from surface 

193 reflectances from the Multi-Angle Implementation of Atmospheric Correction (MAIAC) 

194 algorithm (Lyapustin et al., 201 la, b). We focus on the Aqua GIMMS NDVI in the main 

195 text and show comparable results with the other NDVI data sets in the appendix. The 

196 Aqua satellite has an ascending node equator crossing near 13:30 LT. We estimate errors 

197 as sum of the standard error of the mean and a nominal empirically estimated constant 

198 error of 0.03. 

199 

200 2.3 MERRA reanalysis data 

201 We examine several meteorological fields from the NASA Global Modeling and 

202 Assimilation Office (GMAO) Goddard Earth Observing System Data Assimilation 

203 System version 5 (GEOS-5) Modern-Era Retrospective Analysis for Research and 

204 Applications (MERRA) data set (Rienecker et al., 2011). These include surface skin 

205 temperatures (T s ki n ) and total profile soil wetness (soil moisture), which are from the 

206 Incremental Analysis Updates 2D simulated land surface diagnostics product. We also 

207 use temperature at 2 m above the displacement height and 2 m specific humidity from the 

208 IAU 2D atmospheric single-level diagnostics product to calculate vapor pressure deficit 

209 (VPD, the difference between the actual and saturation-vapor pressure). Here, we use 

210 daily-averaged fields generated at 2/3° longitude by 1/2° latitude resolution. Near-surface 
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211 specific humidity anomalies (Fig. 2a) in July show significantly drier than average 

212 conditions (13-32%) over a large part of the area examined. The T s ki n anomalies show 

213 that the heat wave (up to ~7K above normal for the monthly average) was confined to a 

214 smaller area in the western part of the region (Fig. 2b). Figure 2c indicates that VPD 

215 anomalies are heavily controlled by temperature; the VPD anomalies in August are 

216 smaller than those in July. Soil moisture for both months shows negative anomalies for 

217 all six boxes (Fig. 2d). 

218 

219 2.4 Catchment-CN land surface model simulations 

220 We examine several variables obtained from an off-line run of the Catchment-CN 

221 land surface model (Koster et al., 2014). The Catchment-CN land surface model is 

222 essentially a merger of the energy and water budget framework of the NASA Global 

223 Modeling and Assimilation Office’s Catchment model (Koster et al., 2000) with the 

224 prognostic carbon elements (and thus prognostic phenology elements) of the National 

225 Center for Atmospheric Research/Department of Energy (NCAR/DOE) Community Land 

226 Model 4 (CLM4) dynamic vegetation model (Thornton et al., 2009; Oleson et al., 2010). 

227 The merged Catchment-CN model has some unique features, including the ability to 

228 represent multiple vegetation regimes within a surface element, each static vegetation 

229 regime associated with a different dynamic hydrological regime. The fractional areas 

230 occupied by individual plant functional types in the merged system do not change, but 

231 vegetation growth, soil hetero trophic activity, carbon stocks, and other ecosystem states 

232 (such as those associated with leaf area index) vary prognostically. Comparison of 

233 simulated fPAR with satellite-based estimates from the GIMMS AVHRR dataset (Tucker 
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234 et al., 2005) demonstrate that the model, while biased, captures well the controls imposed 

235 by water supply on the global distributions of phenological variables (Koster et al. 2014); 

236 overall, Catchment-CN is found to be a useful tool for the analysis of the connections 

237 between climate and vegetation. 

238 Fluorescence was added to the model by including the approach detailed for a 

239 similar implementation within the CLM4 (Lee et al, 2014). The fluorescence code uses 

240 as inputs the photosynthesis rate, the intracellular leaf CO 2 , and the CO 2 compensation 

241 point; it produces as an output SIF, as a daily mean for both the sunlit and shaded 

242 portions of the canopy. We used a model calibrated to leaf scale measurements of 

243 chlorophyll fluorescence from pulse amplitude modulated (PAM) fluorometry to simulate 

244 0f as a function of the rate of photosynthesis simulated within the model. Key model 

245 variables are the flux of absorbed PAR, the rate of photosynthetic electron transport 

246 provided by the photosynthesis parameterization, and the level of non-photochemical 

247 quenching that can be measured with PAM fluorometry. 

248 The offline Catchment-CN simulations are driven with atmospheric forcing from 

249 the MERRA-Land reanalysis product (Reichle et al., 2011), which is identical to that of 

250 MERRA except that surface precipitation is corrected to a global, daily, 0.5° gauge 

251 product. Full Catchment-CN model spin-up was ensured by cycling over a 35 year 

252 period several times prior to producing the simulation data examined here. The 

253 Catchment-CN model was run on 64,770 irregularly shaped tiles (or computational 

254 elements) based on watershed delineations with a mean area of 2,010 km 2 and median 

255 area of 1,186 km 2 . We use monthly mean output generated at 2.5° longitude by 2.0° 

256 latitude resolution from 2007 to 2013 for all parameters examined including surface skin 
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257 temperature, fPAR (calculated as APAR/PAR), PAR, SIF, and GPP. 

258 

259 3. Results 

260 3.1 Seasonal anomalies 

261 The top half of each of the six panels of Figure 3 (one panel for each of the box 

262 regions in Fig. 1) shows the climatological seasonal cycles of GOME-2 SIF and GIMMS 

263 NDVI as well as the values of 2010 for April to September. The bottom half of each 

264 panel shows the 2010 anomalies of VPD and soil moisture from MERRA. We next 

265 discuss results for boxes grouped by dominant vegetation types. 

266 

267 3.1.1. Croplands 

268 For the boxes dominated by croplands (1,2, and 5), climatological SIF and NDVI 

269 reach their maxima in middle June to late July depending on location; croplands towards 

270 the east generally peak later. As has been shown in other studies for croplands (as well as 

271 mixed forest), SIF starts to decline earlier in autumn as compared with reflectance-based 

272 indices such as the NDVI; the earlier decline of SIF is in better agreement with GPP from 

273 flux tower measurements (Joiner et al., 2014). Soil moisture anomalies indicate 

274 substantially drier than normal conditions starting around the middle of May for these 

275 boxes. VPD anomalies are large for box 1 that is within the area impacted by the heat 

276 wave. Similar to the surface temperatures, VPD anomalies peak in late July. The SIF and 

277 NDVI 2010 negative anomalies in these boxes are significant. For box 1, within the heat 

278 wave region, there is a somewhat later and smaller 2010 anomaly as compared with the 

279 other two cropland-dominated boxes. This could be because box 1 is in the basin of the 

280 Volga river that supplies ground water. Both SIF and NDVI indicate a slight partial 
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281 recovery in August in boxes 1 and 2 only. There is a very strong correspondence between 

282 the GIMMS NDVI and SIF for all areas. 

283 

284 3.1.2. Grasslands and mixed forest 

285 Boxes 3, 4, and 6 are primarily covered by grasslands and mixed forest. Box 3, 

286 composed primarily of mixed forest, appears to be less affected by drought than the other 

287 regions examined; the differences between the climatology and 2010 for SIF and NDVI 

288 are not statistically significant. Box 4, which is primarily grasslands, shows negative 

289 2010 anomalies for both NDVI and SIF starting in early June. In contrast, box 6, which 

290 contains a mixture of grasslands and mixed forests, shows only small negative 2010 

291 anomalies starting in late June. 

292 

293 3.2 Land surface modeling results 

294 Figure 4 shows monthly means of the Catchment-CN land surface model output 

295 for the climatology and for 2010 in the six boxes. Parameters examined are fPAR, PAR 

296 APAR, and LUE. In Figure 5, SIF and GPP as well as SIF and GPP normalized with 

297 respect to incoming PAR are shown. Surface skin temperature and soil moisture (root 

298 zone) are shown in Figure 6. There are significant negative 2010 anomalies in GPP for all 

299 boxes starting mostly in June, which are influenced by negative anomalies in LUE. The 

300 surface skin temperatures are generally higher in 2010 for all regions as may be expected 

301 in conjunction with lower GPP. Soil moisture shows clear negative 2010 anomalies after 

302 May-April in the most of boxes except boxes 1 and 3, where there are negative anomalies 

303 for the 2010 growing season. 

304 In contrast, the model’s fPAR does not show a 2010 anomaly for box 3 dominated 
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305 by mixed forest. In addition, the model’s fPAR negative anomalies for the other boxes 

306 generally begin in July or August, about one month later than the GPP anomalies. PAR 

307 2010 anomalies, on the other hand, are generally insignificant to positive, owing to 

308 decreases in cloudiness during the peak drought months. Because the model’s 2010 PAR 

309 and fPAR anomalies are of opposite sign, this leads to smaller negative or insignificant 

310 2010 APAR anomalies as compared with fPAR anomalies. 

311 The model’s 2010 SIF anomalies are somewhat smaller (in a percentage sense) 

312 than those of GPP. For example, when GPP drops to near zero starting in August for box 

313 5, while the simulated SIF remains slightly above zero for August 2010. The model’s 

314 2010 SIF anomalies in most boxes are significant starting in July, while GPP negative 

315 anomalies begin in June for all boxes. However, when nonnalized with respect to 

316 incoming PAR, SIF shows earlier negative anomalies (starting in June) for most boxes 

317 and significant anomalies for all boxes, which is similar to the GOME-2 SIF anomalies. 

318 However, the model’s 2010 negative fPAR anomalies start later (July), while the GIMMS 

319 NDVI anomalies begin earlier similar to the SIF anomalies. This indicates that the 

320 model’s fPAR response to drought/heat stress may have occurred somewhat late. 

321 In our analysis of GOME-2 SIF in Fig. 3, we partially filtered for clouds; we 

322 removed pixels with effective cloud fractions > 0. 15. We also nonnalize SIF with respect 

323 to the incoming clear-sky PAR. It should be noted that the spectral signature of SIF is not 

324 affected by clouds. The main effect of clouds on satellite-observed SIF is a shielding 

325 effect that reduces the amount of canopy-level SIF that is observed by the satellite 

326 instrument. The cloud-shielding effect is relatively small for thin and broken clouds with 

327 low cloud fractions. For example, Frankenberg et al. (2013) showed with simulated data 
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328 that 20% or less of the canopy-level SIF signal is lost from satellite observation for cloud 

329 optical thicknesses up to 5. To be consistent, because the PAR-normalized, cloud-filtered 

330 GOME -2 SIF is biased toward clear skies, it should be compared with the PAR- 

331 nonnalized SIF from the model. 

332 The model produces similar (PAR-normalized) SIF anomalies as compared with 

333 the GOME-2 data, although the overall phenology is somewhat different. One difference 

334 between model and GOME-2 SIF is for the mixed forest dominated box 3. For this box, 

335 GOME -2 SIF does not show a significant 2010 SIF negative anomaly, while the model 

336 simulates a significant (nonnalized) anomaly. The fact that NDVI does not show a 

337 significant 2010 anomaly for this box is consistent with the absence of an fPAR anomaly 

338 in the model. Therefore, the model’s negative 2010 photosynthesis anomaly may be 

339 overestimated for this box. 

340 To provide an overall regional context, Figure 7 shows maps of 2010 anomalies 

341 of fPAR, PAR, APAR, and LUE from the land-surface model for July and August. fPAR 

342 anomalies are smaller in July as compared with August. The higher positive 2010 PAR 

343 anomalies in July are reflected in the APAR anomalies and lead to some positive 

344 anomalies in APAR. LUE anomalies are negative over most of the domain and more 

345 significant in August. 

346 Figure 8 also shows maps of 2010 anomalies from the land-surface model for SIF, 

347 GPP, and both quantities normalized with respect to PAR. As noted above, the positive 

348 anomalies in GPP and SIF in the northwestern portion of the study area result from PAR 

349 anomalies, while the negative anomalies towards the south in the PAR-adjusted quantities 

350 are shown with contributions from fPAR; the increase in magnitude of the negative 
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351 anomalies from July to August results primarily from the decline in fPAR over that 

352 period. 

353 Figure 9 compares the model’s SIF with that from GOME-2 for both the 

354 climatology and 2010 anomalies. Here, the model SIF is nonnalized with respect to PAR 

355 and GOME-2 SIF is scaled as before by cosine of the solar zenith angle. The satellite SIF 

356 data are shown at both the model resolution and a higher spatial resolution. To provide 

357 more samples per gridbox, we retain all data with effective cloud fractions up to 0.3. This 

358 did not substantially change the spatial or temporal SIF distributions as compared with a 

359 lower cloud fraction threshold. The satellite and model SIF (climatology and anomalies) 

360 are generally comparable, although there are some differences in the spatio-temporal 

361 distributions. Overall, the model is shown to produce a reasonable response of SIF to the 

362 drought/heat wave. At the same time, it provides insight into how the different 

363 components of SIF and SIF itself may respond to heat and water stress. Note that the 

364 model data are output as monthly means (averages of daily means) and so cannot be 

365 directly compared with instantaneous satellite SIF measurements taken at a specific time 

366 of day. 

367 

368 3.3 Inter-annual variations in SIF and NDVI 

369 Figure 10 compares interannual variability (2007-2013) of the GOME-2 SIF and 

370 the GIMMS NDVI integrated over April-September for the six boxes examined above. 

371 Note that the axes are normalized to the maximum values. For all boxes except box 3, 

372 SIF and NDVI are correlated (r 2 values of 0.75-0.91). This relatively high correlation 

373 confirms that fPAR is a major contributor to the interannual variability of SIF in this 

374 region. 
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375 An interesting feature is the deviation of the fitted slopes (solid lines) from the 

376 one-to-one (1:1) lines (dashed). For example, for box 4 (primarily grasslands), the 

377 minimum value of SIF in 2010 is > 60% less than the maximum, while that of NDVI is 

378 -35% less than maximum. While fPAR impacts both SIF and NDVI, SIF is additionally 

379 affected by fluorescence efficiency, related to photosynthesis and light-use efficiency. 

380 This may explain the larger percentage drought impact on SIF as compared with NDVI 

381 for these boxes. It should also be noted that fPAR is somewhat non-linear with respect to 

382 NDVI ( e.g Los et al., 2000). 

383 

384 4. Conclusions 

385 We have examined the response of canopy-level SIF to heat and drought stress in 

386 2010 over a portion of Russia that includes both agricultural areas and forested regions 

387 using satellite SIF and NDVI observations as well as model simulations. SIF and NDVI 

388 satellite data show similar signs of drought stress early in the growing season well before 

389 the onset of the heat wave both inside and outside the main area of the heat wave. Large 

390 declines in 2010 are seen in both quantities throughout much of the drought-affected area. 

391 Areas dominated by crops and grasslands showed significant drops in SIF and NDVI, 

392 while regions of predominantly mixed forest showed small to insignificant reductions. 

393 We simulated SIF using a global land surface model forced by observations-based 

394 meteorological fields. The model simulated large negative anomalies in 2010 SIF similar 

395 to those seen in the GOME-2 satellite SIF data. The model also produced spatial and 

396 temporal patterns of the SIF anomalies similar to those derived from GOME-2, although 

397 with some exceptions. There exists potential to improve the model’s response by using 
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398 the satellite SIF observations for data assimilation (modification of the model’s 

399 prognostic variables) and/or parameter estimation; this could be a topic of a future study. 

400 Although the model simulated earlier drought-related declines in photosynthesis as 

401 compared with fPAR, the NDVI data suggest that there were significant declines in fPAR 

402 early in the growing season for areas dominated by crops and grasslands. 

403 New satellite sensors, such as the recently launched Orbiting Carbon Observatory 

404 2 (OCO-2) (Frankenberg et al., 2014) and the TROPOspheric Monitoring Instrument 

405 (TROPOMI) (Veefkind et al., 2012) to be launched in 2016 will offer higher spatial 

406 resolution measurements as compared with GOME-2. In addition, these satellites will 

407 make measurements from sun-synchronous polar orbits with local overpass times in the 

408 early afternoon, when stress effects should be peaking and may be larger during the 

409 morning overpass of GOME-2. We plan to utilize these new data sets for further 

410 examination of the manifestation of stress effects on observed SIF. We also plan further 

411 comparisons between satellite and modeled SIF with an aim towards using the satellite 

412 SIF data to improve models as demonstrated by the pioneering study of Zhang et al. 

413 (2014). 

414 

415 Appendix 

416 Here, we compare seasonal cycle of NDVI from GIMMS, MAIAC, and the 

417 standard Aqua MODIS product MYD13Q1, 16-day L3 global 250m SIN grid collection 5 

418 (MYD13) for the climatology and 2010. The products differ mainly in how the cloud 

419 detection is applied. The MYD13 data have been additionally filtered for cloud and 

420 aerosol contamination following the methodology of Xu et al. (2011). All three NDVI 
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421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 


data sets look similar, although there are a few exceptions. For example, in box 1, the 
climatologies from GIMMS and MAI AC are similar, but MYD13 shows a significantly 
lower peak than the other two. Also, climatological MYD13 in box 4 does not show a dip 
in middle June as the GIMMS and MAIAC do. For other boxes, the seasonal cycles for 
all three products are more similar. 


GIMMS MAIAC MYD13 







Figure Al: Seasonal cycles of NDVI from GIMMS (blue), MAIAC (green) and standard product 
of MYD13Q1 (red); solid lines (broken lines with symbols) are for the climatology, (2010). 
Averages are computed using data only where all three data sets provided successful retrievals. 
MYD13 data are interpolated to match 8-day intervals of the other data sets. 
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661 Table 1 : Location of the twelve box regions (center of 2° longitude x 1° latitude box) and 

662 major IGBP vegetation types; CRO: croplands, GRA: grasslands, MF: mixed forest. The 

663 percentage of the coverage is also shown (not shown if the coverage < 5%). 


Box number 

Latitude 

Longitude 

Vegetation cover (%) 

1 

54.5° N 

48.0° E 

CRO: 62 
GRA + MF: 14 
MF: 11 

2 

52.5° N 

55.0° E 

CRO: 97 

3 

54.5° N 

57.7° E 

MF: 95 
GRA + MF: 5 

4 

51.0° N 

67.5° E 

GRA: 100 

5 

54.0° N 

67.5° E 

CRO: 81 
GRA: 10 
GRA + MF: 8 

6 

56.5° N 

70.5° E 

GRA + MF: 86 
MF: 8 
CRO: 4 
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666 Figure 1 : Map of land cover type for 2010. The six box regions used for further analysis 

667 are also shown. 
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Figure 2: Maps of July (left column) and August (right) 2010 anomalies of MERRA 
meteorological fields (differences between July (August) 2010 and average of all other 
July’s (August’s) from 2007-2013 not including 2010): a) specific humidity (anomalies 
in terms of %), b) surface skin temperature (anomalies in K), and c) vapor pressure 
deficit (anomalies in terms of kPa). 
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Figure 3: Top panels: Seasonal cycle (8-day means) of GOME -2 fluorescence 
[mW/m 2 /nm/sr] (black lines), Aqua MODIS GIMMS NDVI [unitless] (magenta lines); 
solid lines (broken lines with symbols) are for climatology (2010). Error ranges are 
indicated as shading or vertical bars where for clarity only a few representative error bars 
are shown (in July) for the 2010 data. Bottom panels: Vapor pressure deficit anomaly 
[hPa] (black line) and soil moisture (SM) anomaly [fraction] (blue line) for the six boxes 
shown in Fig. 1. Anomalies are calculated as 2010 - climatology. 
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Figure 4: Monthly mean Catchment-CN land surface model results with MERRA forcing 
for selected boxes shown in Fig. 1. Black (red) lines represent climatological mean values 
(2010 values). From left column, fPAR [unitless], PAR [W/m 2 ], APAR [W/m 2 ] and d) 
LUE [ug C/J]. LUE is calculated as GPP/ APAR. The black vertical bars indicate standard 
deviations. The blue vertical dotted lines indicate June. Note: different y-scales are used 
for the different boxes. 
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Figure 5: The same as Fig. 4 but for (from left) SIF [nmol photons/m 2 /s], SIF normalized 
with respect to PAR [x 10" 3 ], GPP [jig C/m 2 /s] and GPP normalized with respect to PAR 
[p,g C-W/s]. 
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Figure 6: Same as Fig. 4 but for surface skin temperature [K] (left) and root-zone soil 
moisture [fraction] (right). 
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Figure 7: Maps of 2010 anomalies for July (right column) and August (left), computed as 
differences between July (August) 2010 and average of all other July’s (August’s) from 
2007-2013 not including 2010 calculated using MERRA-forced land surface model 
simulations: a) fPAR [unitless], b) PAR [W/m 2 ], c) APAR [W/m 2 ] and d) LUE [ug C/J]. 
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Figure 8: Same as Fig. 7 but for: a) SIF [|imol photons/m 2 /s], b) SIF normalized with 
respect to PAR [x 10" 3 ], c) GPP [ug C/m 2 /s], and d) GPP normalized with respect to PAR 
[p,g C/s/W]. 
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Figure 9: Maps of the SIF monthly climatology (a, b and c) and anomaly (d, e and f) for 
July (left column) and August (right) from MERRA-forced land surface model 
simulations (a and d), GOME-2 with 2.0°x2.5° resolutions (b and e), and GOME-2 with 
0.5°x0.5° resolutions (c and f). Anomalies are computed as in Fig. 3. Model SIF is 
normalized with respect to model PAR. 
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722 Figure 10: Scatter diagram of April-September integrated GOME-2 SIF and Aqua 

723 MODIS GIMMS NDVI for each year in the range 2007 to 2013 for the six boxes shown 

724 in Fig. 1; red numbers indicate years (i.e., 07=2007). Values are scaled (divided by the 

725 maximum for each box). Solid line: linear fit; dashed: 1 : 1 line. The dominant vegetation 

726 type, correlation (r 2 ), and slope values are provided for each box. 

727 

728 
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